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^D , We discuss the production of forward jets in high energy processes 

where one probes a dense hadronic wavefunction. In particular, and as 
r^ , a signature of parton saturation, we discuss the possibility of a strong 

pi i' momentum decorrelation in Mueller- Navelet jets which leads to a geometric 

JL I scaling behavior. 

(D 
^; PACS numbers: ll.15.Kc, 12.38.Cy 

1. Introduction 
> , 

On , In high energy collisions of two hadrons, a large number of the produced 

particles emerge in the forward directions. If such a particle of energy E 
forms an angle 9 with the beam axis, then it is convenient to decompose its 
3-momentum in terms of its rapidity, representing the longitudinal momen- 

f^ ' tum, and its transverse momentum. For massless particles the rapidity is 

equal to the pseudorapidity t] = — lntan(6'/2), while the magnitude of the 
transverse momentum is k± = EsinO. Since the forward directions corre- 
spond to 9 being close to or vr, the standard belief until a few years ago 
was that only soft physics is involved and one would have to resort on non- 
;J] ■ perturbative phenomenological descriptions. However, it is obvious that a 

5r ! large energy of the produced particle can compensate for the smallness of 

9, thus rendering k± a perturbative scale. For instance, in p-p collisions at 
the LHC it could be possible to measure jets at forward rapidities up to 
\r]\ ~ 6.5 and transverse momenta as low as k± ~ 20 GeV. 

In such a kinematic regime it is clear that one can use weak coupling 
techniques. But the largeness of the rapidity indicates that fixed order 
perturbation theory might fail and calls for possible resummations to obtain 
more reliable results. Indeed, this is what sometimes is called the semihard 
region of QCD, where k± is much larger than Aqcd but at the same time 
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much smaller than the total energy ^/s of the process and large logarithms 
of the type In s/k\ appear and need to be resummed. 

Still, the aforementioned resummation might not be enough since when 
looking at forward directions one is sensitive to the the softer, small-x, com- 
ponents of the wavefunction of the one of the incoming hadrons. Here by 
components we refer to quarks and gluons, with the latter dominating be- 
cause of their vectorial nature and the 3-gluon vertex in QCD. In the high 
energy limit and/or for large nuclei (due to an ~ A^'^ enhancement coming 
from the 2D density of nucleons), the hadronic wavefunction becomes dense, 
partons with sufficiently low, but perturbative, k± can overlap and interact 
with each other leading to a state of saturation [T] . Technically this means 
that one has to perform a resummation of logarithms in the presence of a 
background field which in general can be strong. This leads to nonlinear 
equations for the evolution of the partonic densities, with the non- linear 
terms clearly corresponding to this partonic overlap and becoming impor- 
tant in the saturation regime. Pictorially this is shown in Fig. [1] which is 
our standard picture for parton evolution in QCD. DGLAP is the evolution 
corresponding to an increase in transverse momentum, and even though it 
can lead to an increase in the number of partons, these become smaller 
and smaller due to the uncertainty relation, and the wavefunction remains 
dilute. BFKL is the evolution corresponding to an increase in ln(l/x), the 
number of partons again increases but now they are typically of the same 
size and the wavefunction becomes denser. For a given A;_|_ we reach a point 
in X where gluons start to feel the presence of each other, and these par- 
ticular modes saturate reaching their maximal allowed value of order l/a^. 
Clearly modes with a larger k± will saturate at a smaller value of x, since 
the partons will be smaller in size. Therefore one arrives at the concept 
of the saturation momentum, which increases with ln(l/x), and which is 
the borderline between the saturated and the non-saturated modes of the 
hadronic wavefunction. Some estimates that one can give for the value of 
the saturation momentum are Q^ ~ 1 GeV^ for protons at HERA, Qs ~ 1-5 
GeV^ for gold nuclei at RHIC, Q^ ~ 2 GeV^ for protons at LHC and Q^ ~ 6 
GeV^ for lead nuclei at LHC. 



2. Dijet production at lowest order 

At lowest order in perturbation theory the production of two jets is 
shown in Fig[2j To the order of accuracy the incoming partons, here gluons, 
participating in the process have zero transverse momentum pi^ = p2±_ = 0. 
Momentum conservation along this transverse plane implies that fcij_ + 
fe2± = 0, that is, the azimuthal angle between two outgoing jets is Ac/) = vr. 
Longitudinal momentum and energy conservation constrain the longitudinal 
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Fig. 1. Parton evolution in QCD. 



Fig. 2. Two jets at lowest order. 



momentum fractions xi and X2 of the incoming partons to be expressed in 
terms of the center of mass energy ^/s and the outgoing jets momenta as 



Xl 



-i± 



,Vi 



+ 



hi. 



:.V2 



X2 



^1± 



-VI 



+ 



k2i. 



-m 



Denoting by PS the phase space d^fci_Ld^fci_Ld?7id?72 we have 



(1) 



da 
dPS 



"^ Xlfi{xi,fl^) X2fj[x2,^?) (5^^^ (fcl± + k21_) -77^ 
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where k± is the common magnitude of the momenta, i,j refer to the various 
parton species, xfi{x,iJ,'^) are the corresponding distribution functions and 
(Tjj the partonic cross sections. The emission of an unobserved extra gluon 
or quark will smear a bit the peak, but not much, since it is of higher order 
in perturbation theory and we assume the coupling to be weak. 



3. Two jets in the forward direction 

Because of the forward kinematics fixed order perturbation theory is 
not sufficient. Let us first consider the case where both jets are produced in 
the forward direction as shown in Fig. [3l that is we assume that rji and r/2 
are large, say positive, and close to each other. We have in mind collisions 
of asymmetric objects, like d-Au at RHIC [2], with the forward direction 
being the one of the colliding deuteron. From the kinematics in Eq. ([1]) it 
is straightforward to see that for fixed transverse momenta xi can be close 
to 1, but X2 is necessarily very small. The nuclear wavefunction can be 
dense due to the A^'"^ enhancement and due to the resummation of large 
[ois ln(l/2;2)]"' terms, with a^ = agNc/TT and Nc the number of colors. High 
energy evolution proceeds via gluon cascades, as shown in Fig. [3l which 
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Fig. 3. Pair of jets in forward direction. 



Fig. 4. Mueller-Navelet jets. 



are ordered in longitudinal momentum. Such cascades carry significant 
transverse momenta, during the scattering they are released and thus there 
is a large imbalance in the transverse momenta of the two forward jets and 
the "away peak" at A(j) = vr can disappear. In the BFKL regime one has 



da 
dPS 



^Xi/j(xi,^^)$g(x2,fcl± + k2±) 






(3) 



where $g is the unintegrated gluon distribution function whose presence is 
necessary in order to describe the transverse momentum "built-up" along 
the cascade. Note that Eq. ^ is valid only so long as the gluon density is 
not too high and the non-linearities can be neglected. In the presence of 
saturation the cross-section involves also high point gluon correlators. 



4. Mueller-Navelet jets 

Finally we consider the forward-backward case in p-p (or p-p) collisions, 
that is when the two jets are separated by a large rapidity interval. These are 
the Mueller-Navelet jets [3], where more precisely one has Y = rji — rj2 ^ 1 
with \r]i\ ~ \r]2\ and the corresponding diagram is shown in Fig. [H We 
have decided to view the process in a different way than the one of the 
forward dijet in the previous section. Now the resummation is performed 
in the partonic cross section, since the emitted gluons in the final state are 
ordered in rapidity in order to get the dominant contribution. Nevertheless, 
one can show that the kinematics in Eq. ([1]) holds approximately and we 
have xi^/s ~ ki±e'^^ and X2\/s ~ k2j_e~'^^ , so that one can choose xi and 
X2 to be "large", say around 0.1. We can write the cross section as 
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As explained earlier, at the Born level we get two back-to-back jets and the 
partonic cross section is dij ~ 0{1). An extra "minijet" introduces some 
decorrelation in the transverse plane and iTjj ~ ©(a^^), with two minijets 
there is more decorrelation and dij ~ 0{a^Y'^) and so on. We integrate 
over the minijets phase space and sum for large Y to find (for any i,j) 

^^ - "' ^(l^,fci±,fc2±), (5) 



where $ satisfies the BFKL equation. It grows exponential with Y for fixed 
kix_ and fci_|_ and imposes transverse momentum decorrelations in the angle 
and the magnitude [il lSJH 

As a particular example of these magnitude decorrelations, let us recall 
that BFKL dynamics with the necessary inclusion of unitarity corrections 
leads to geometric scaling [THS], which was introduced and had its greater 
success in e~p DIS at HERA [lOl. Integrating the transverse momenta 
above the thresholds Qi and Q2 we have 

~ ^eff 7^2 / TT^ T^ ^ ^^v[asX{l)Y\ + • • • , (6) 



dxidx2 " Q\ Ji-ioo 27ri 7(1 - 7) 

where -Fefi stands for the contributions from all the quark, antiquark and 
gluon collinear distributions of the two colliding protons together with some 
color factors, x{l) is the known eigenvalue function of the BFKL equation 
and the dots stand for unitarity corrections. When Q\ is much larger Q2^ 
but in such a way that the integrand in Eq. ([6]) varies slowly, the integration 
in the saddle point approximations leads to [11] 



d^ _p. ^^QiV^ .7^ 

In the above Ql = Q^i^^^^ ~^'''> is a "saturation scale", with A^ ~ 0.3 [T2], 
Yq ~ (l/os) ln(l/as) and 7^ = 0.372 is the anomalous dimension dominat- 
ing the integration. Eq. d?]) exhibits geometric scaling since, leaving aside 
the slowly varying prefactor Fcg, it depends only on two variables, Q\ and 
the combined one Q\IQ\-, instead of the three Qf, Q\ and Y . This is analo- 
gous to the scaling observed in 7*p DIS, with the analogy being more precise 
if we let Q\ — )• I^q.qy> ^^d Qi ^ Q^, with Q^ the photon virtuality. Thus the 
softer jet looks like the "target" with the harder jet being the "projectile". 
The difference with DIS is that here we have a very large saturation scale 
to start with, which is by definition the threshold transverse momentum of 



^ An NLO BFKL calculation for some observables leads to results that are very close 
to those obtained from an NLO DGLAP analysis |6]. 
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the softer jet and which furthermore evolves to higher values with increasing 
Y. The price to pay is that the cross section is proportional to the "target" 
size squared which is l/Qi- The rapidity Yq corresponds to the amount of 
evolution which is needed to saturate a small hadronic system, like a dipole 
or a high-momentum parton, on the resolution scale set by its own size. An 
NLO estimate gives Yq — 8 which is very well within the reach of the LHC. 
According to Eq. d?]), one can explore the scaling behavior of this dijet 
cross section by keeping the transverse and longitudinal momenta of the 
softer jet fixed, and vary those of the harder jet so that xi remains fixed. 
As said, this represents a strong momentum decorrelation and would test 
BFKL dynamics and saturation. 
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